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Abstract 

The accurate prediction of the neutrino beam produced in muon decays and 
the absence of opposite helicity contamination for a particular neutrino flavor 
make a future neutrino factory the ideal place to look for the lepton flavor vi- 
olating (LFV) decays of the kind fi + — > e + v e v fJi and lepton number violating 
(LNV) processes like fi~ — ► e~v e y il . Excellent sensitivities can be achieved using 
a detector capable of muon and/or electron identification with charge discrimi- 
nation. This would allow to set experimental limits that improve current ones 
by more than two orders of magnitude and test the hypothesis that the LSND 
excess is due to such anomalous decays, rather than neutrino flavor oscillations 
in vacuum. 
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1 Introduction 



The Standard electroweak Model (SM) is built on an absolute conservation of 
three separate lepton flavors (LF): the electron flavor L e , the muon flavor L„, 
and the tau flavor L T . The total lepton number is defined as Ltot = L e + + L r . 

In a model independent way, the main decay mode of the positive muon[] into 
a positron and two neutrinos can be written as|]l|, 0] : 

fi + ^e + + n + n' (1) 

where n and n' denote neutrinos, which can be either neutrinos or antineutrinos 
and of any flavor e, fi or r. In the SM, individual LF conservation implies that 
n = u e and n' = v^. 

The evidence from neutrino oscillation searches that the neutrinos are in fact 
massive and mixed, implies that the LF conservation is not exact. However, LF 
conservation is constrained by stringent experimental limits obtained in processes 
involving charged leptons. For example, the present experimental 90%C.L. upper 
limits on the most interesting of these decays are: 

(2) 
(3) 
(4) 
(5) 

These stringent limits are not inconsistent with the neutrino oscillation results 
since given reasonable upper bounds on the neutrino masses, the effect induced 
within the SM would be too small to be seen]?]]. However, within extensions 
of the SM, neutrino oscillations raised the possible prospect that there might 
exist observable processes that violate the charged-lepton number p. Projects 
are currently underway to improve several of these upper limits significantly. 

In contrast to this, direct experimental limits on LF conservation in decays 
involving neutrinos are much less stringent. For example, the limit 

BR(p~ -» e~v e Vp) < 1.2 x 1(T 2 : [9] (6) 

is orders of magnitude worse than limits involving charged leptons. In extensions 
of the SM, the existence of LF violation could open new modes of decay where 
n and n' in Eq. (jl[) are indeed any given neutrinos or antineutrinos. 

One would clearly like to experimentally test these limits by a sensitive study 
of the neutrinos produced in such decays. Such studies could directly be relevant 
to a further understanding of the neutrino sector, as described in the following 
sections and would represent complementary investigations to those for neutrino 
flavor oscillations. 
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If the excess of events found in the Los Alamos LSND experiment [10, 11] is 



interpreted as due to fi + — > e + + v e + n (with a branching ratio equal to the 



1 Similar arguments hold for negative muons. 
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measured probability of v e appearance in case of flavor neutrino oscillations), 
this anomalous decay could be easily tested at a future neutrino factory with a 
detector capable of charge discrimination, looking at the more convenient decay 
pT — ► e~v e y iii as described in section |7[ 

Good sensitivities could also be reached at a neutrino factory for the search 
/i + — > e^vgVu described in section ^ in which one look for interactions of Un giving 
fi" in a pure beam that produces /i + . 



2 The neutrino oscillation sector 

Hints that neutrinos are massive particles come from the observation of three 
anomalous effects — the LSND excess [|To|, 11], the atmospheric anomaly [13 



14, |l5|, 16, 17] and the solar neutrino deficit [18, IS, ^(| |2l|, |^]. In particular, 
the atmospheric results are the most convincing ones. All three effects can be 
naturally explained in terms of neutrino flavor oscillations, which will occur when 
neutrinos propagate through space, if their masses are non-degenerate and the 



weak and mass eigenstates are mixed [ 23 1. 



However, in order to explain all three experimentally observed effects in terms 
of neutrino flavor oscillations, one is forced to invoke additional sterile neutrino 
states [^4]] to accommodate the very different frequencies of oscillations — given 
by the mass differences squared Am 2 's — indicated by the three different effects. 
The existence of such neutrinos is a currently unresolved problem and clearly 
demonstrates that the neutrino sector is not fully understood. Several attempts 
were made to explain all data in terms of only three massive neutrinos j2l|. But 
they are all excluded by the latest data. 

From a phenomenological point of view, we recall that the neutrino flavor 
oscillation hypothesis predicts a well defined dependence of the phenomenon as 
a function of the neutrino energy, characterized by the so-called L/E behavior, 
where L is the distance between the source and detector and E the neutrino 
energy. So far, no experiment has conclusively demonstrated such a L/E depen- 
dence of the anomalous effect, with maybe the exception of the SuperKamiokande 



data which favors [26] a dependence oc LE n where n ~ — 1. 

In such an unclear situation, is it possible to envisage "non-flavor-oscillation" 
mechanisms to explain part of the neutrino data? 

Aside from theoretical arguments against sterile neutrinos, we argue that, 
from a phenomenological point of view, the LSND effect is particular: it has 
a small probability, measured to be (2.5 ± 0.6 ± 0.4) x 10~ 3 [11], in contrast to 
the solar and atmospheric neutrino anomalies, which are large. Hence, LSND 
is a natural candidate for an interpretation involving a different physics than in 
atmospheric and solar neutrino flavor oscillations. 



3 Relevance to the LSND puzzle 

We recall that the LSND effect was first reported as an excess of iVs in the 



f^'s flux from the /i + Decay- At-Rest (DAR) process] 1C]. The neutrino beam 
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is obtained with 800 MeV kinetic energy protons hitting a series of targets, 
producing secondary pions. Most of the tt + come to rest and decay through the 
sequence ir + — ► n + v u , followed by fi + — » e + i/ e i7 (U , supplying the experiment with 
the P-'s with a maximum energy of 52.8 MeV. The intrinsic contamination of 
P e 's coming from the symmetrical decay chain starting with ir~ is estimated to 
be small since most negatively charged mesons are captured before they decay. 

The excess of P e 's, explained in terms of neutrino flavor transitions of the 
type P/j, — ► P e , occurs via the reactions: 

4. 4- _ _ vacuum _ _ 4- /— \ 

-> e T v e ^; i/ M — > i/ e ; u e p ^ e^n (7) 
Additional evidence in favor of neutrino oscillation was reported in the Decay- 



In-Flight (DIF) sample, though with a low statistical significance [12]. In a recent 
re-analysis of the complete data sample collected [|Tl|], the significance of the DIF 
data seems to be even lower. Hence, we concentrate on the hint from stopped 
muons, and ignore the DIF result. 



The latest KARMEN2 results [£7| come very close to contradicting the LSND 



claim, however the experimental sensitivity is marginal to conclusively exclude or 



confirm completely the LSND excess. A new experiment, MiniBOONE[28], will 
confront the flavor oscillation hypothesis with a very high statistical accuracy. 

In case of a negative result, one will only be able to conclude that the LSND 
excess was not due to neutrino flavor oscillations. 

The implications of exotic muon decays on the LSND excess has been studied |2J 
showing that two explicit models predicted interactions of about one order of 
magnitude smaller than what would be relevant for LSND. In a model indepen- 
dent proof in which any contributions from neutrino mixing is neglected, the 
authors of Ref . [|30|| prove that new lepton flavor violating interactions, under the 
constraint of LF data involving charged leptons, fail short to explain the LSND 
effect by a rate factor of almost three. However in Ref. [31], it is reported that 
exotic decays that produce two antineutrinos 

li + -> e + + P e + Pi (£ = e, n, r) (8) 

cannot be ruled out as the cause of the LSND excess in a model independent 
way. 

Regardless of any theoretical prejudice, the excess of electrons found by the 
LSND experiment stands today as a still unresolved puzzle of neutrino physics. 



4 Searches at a neutrino factory 

A neutrino factory^, [jjj is understood as a machine where low energy muons of 
a given charge are accelerated in a storage ring. The two neutrinos n, n' produced 
in the decay Eq. (|l]) will be boosted in the forward direction of the muon flight 
path. Hence, the muon storage ring is composed of long straight-lines in order 
to produce directional neutrino beams. 

In such machines, muons are produced in decays of secondary pions produced 
by few GeV protons incident on a target. In current designs, muons are captured 
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with high efficiency and very high integrated protons-on-target intensities are 
envisaged in order to produce very intense neutrino sources. 

The neutrino physics potentialities of such machines has been largely dis- 
cussed in the literature |3f|. In particular, we mention our study in the con- 
text of a short-baseline experiment to search for neutrino flavor oscillation in a 



background free environment in the Am 2 region indicated by LSND[36|. 

A neutrino factory is also an ideal place to study neutrinos from exotic muon 
decays. The envisaged flux of neutrinos is sufficiently high to obtain large statis- 
tics of neutrino interaction events. More importantly, we can take advantage 
from the fact that the neutrino beam is produced from muons of a definite sign 
and therefore the decay processes can be studied with a pure initial state. This is 
not the case for traditional pion decay neutrino beams, in which contaminations 
are always present at some level. 

The flavor of the interacting neutrinos can be tested via their charged current 
processes. In case of purely lepton flavor conserving decays /U + — > e + v e i>^, we 
expect to detect only 

v e + N -> e~ + X (9) 
+ N -> /i+ + X (10) 

while exotic decays can be immediately identified by various processes 
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where n stands for neutrinos or anti-neutrinos. 

Charge discrimination of electrons and muons can trivially separate the two 
types of decays. The presence of taus is more difficult to identify but can be 
achieved using topological or kinematical signatures. It requires neutrino beams 
of high energy in order to exceed the tau lepton production threshold. We do 
not consider tau identification any further in this paper, and concentrate on the 
identification of electrons and muons in a low energy setup. 

In order to predict the energy distribution of the (anti)neutrinos in the de- 
tector we assume for definiteness two types of generic decays: 

(i + -> e + + ve + (15) 
/i~ — > e~ + u e + vi (16) 



We address the search for LFV, Eq. (|15[), using wrong sign muons, which are 
experimentally simpler to detect. In order to profit from the enhanced cross- 
section of neutrinos versus antineutrinos, it is better to select positive muons in 
the storage ring, since in this case, the LFV decays produce i/^'s. 

In case we consider the LNV decay of Eq. (|l6|), for which AL = 2, we should 
look for wrong sign electrons. We select yT in the ring since the signal searched 
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for in this case has two neutrinos in the final state, therefore we profit from 
the enhanced neutrino cross sections. Naturally, both signs of muons could be 
studied in a real experiment, in order to provide possible checks for different 
behaviors in n + or \i~ decays. 

Following the discussion in LNV interactions can naturally arise via 
mixing of heavy bosons that transform differently under the Standard Model 
group but identically under the unbroken U(1)em- In particular, the effective 
four-fermion operator relevant for the reaction: 

Ml ->■ e R + u e + Vt (17) 

which violates lepton number conservation by two units, has the form (fJ,LV e )(v£en). 
It couples v e to fi~ and e~ to vg. This operator can be induced, for example in 
supersymmetric models without R-parity, through the mixing of sfermions that 
are SU(2)l singlets with sfermions that are SU(2)l doublets. 

After explicit calculation, we obtain that the square of the scattering am- 
plitude for the decay (|l7|) is < \A4\ 2 > oc (p^ ■ p Ue )(Pe • Pv^) (where pi is the 
four-momentum of particle i). This expression is similar to the one obtained for 
the standard decay n~ — ► e~P e f M ; therefore the flux of f e 's coming from ( |l7|) is 
equal to the one of P e 's produced in the standard fj,~ decay as shown in Figure ||. 

5 Experimental considerations 

We address a few experimental considerations in the context of an optimization 
of searches for LFV and LNV decays of the muon. 

5.1 Beam setup 

Unlike for neutrino flavor oscillations in vacuum, the distance L between source 
and detector is in this case an irrelevant physical parameter. It is hence advanta- 
geous to place the experiment close to the source in order to gain flux like l/L 2 
due to the beam divergence. 

We think that a neutrino beam energy of a few GeV would be optimal in 
order to facilitate the discrimination of the muon and electron charges, and in 
order to reduce misidentified electrons or muons background coming mostly from 
neutral current interactions. We therefore consider a low-energy muon storage 
ringp] with three possible muon energies = 1, 2 or 5 GeV. 

In Table p], we hst the expected event rates from standard muon decays per 
ton of target and 10 19 standard muon decays. The detector is located at a 
distance L = 100 m from a 100 m long straight section of the storage ring. The 
expected event energy spectra are shown in Figure 

2 We also note that studies indicate that the cost of a neutrino factory is driven by the muon energy 
and rises very rapidly with E^, hence a low energy muon beam is also financially favored. 
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Table 1: Expected event rates (CC=charged current, NC=neutral currents) per ton of 
target per 10 19 standard muon decays in a storage ring with a straight line of 100 m and 
located 100 m away from the neutrino detector. Only 50% of the muons are assumed 
to decay in the direction of the detector, the other 50% are lost. 



5.2 Detector parameters 

Experimentally, the presence of LFV decays will characterize themselves by the 
observation of events with "wrong sign leptons". It is therefore mandatory to 
ensure a very good and efficient determination of the lepton charge. 

The muon charge is most easily determined with the help of bending in a 
magnetic field. The radius of curvature in meters in a 1 T field is approximately 
10p(GeV)/3, or 3 meters for p = 1 GeV. In the case of electrons, the capability 
to measure the charge is limited by the radiation length of the target which 
determines the distance after which the electromagnetic shower develops at a 
level where the primary electron is not distinguishable any longer. 

At low energies, we expect events to exhibit simple topologies and many of 
them will be quasi-elastic-like. Hence, the events will be dominated by a leading 
hard lepton accompanied by few soft hadrons. Such events were traditionally 
best studied in bubble chambers, due to the required low density, high granularity 
and homogeneity to capture soft escaping tracks at all angles. In addition, a low 
density, high granularity instrumented target is mandatory to efficiently recognize 
electrons and to discriminate muons from pions. Different detector configurations 
meeting the previously mentioned specifications can be clearly envisaged. 



6 Results for —> e + + U£ + 

We consider a 10 ton fiducial mass detector located at a distance of 100 m from 
the muon storage ring. 

For a 2 GeV muon ring energy, the expected event samples for a total of 10 19 
standard /i + decays are 72'000 CC and 215'000 u e CC events, and 18'300 
NC and 36'600 v e NC. 

To estimate the signal efficiency, we assume that the LFV decay proceeds 
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Table 2: Effect of cuts on background and signal. We assumed a positive muon ring 
energy E^+ of 1, 2 and 5 GeV and a total of 10 19 standard decays. The LFV decay 
has the branching probability of 2.5 x 10~ 3 . Backgrounds come from hadrons escaping 
the detector without interacting or muons from meson decays. For the charged current 
background, no veto on the positive muon has been included. E^- is the energy of the 
identified negative muon in the event. 



through a similar diagram as the standard muon decay as given by the V — A 
theory, however, with interchanged neutrino flavors. We can then essentially 
assume that at the detector location the flux of v^s from LFV decays is similar 
to the one of z/ e 's in the standard muon decay. It is clear that other type of 
interactions could be envisaged and could lead to different energy spectra which 
can be experimentally tested by studying the visible energy distribution of wrong 
sign muon events. 

If LFV decay occur with a branching Br (LFV), we expect to observe the 
number of negative muon events 

N^-,cc = f <S>LFv(v^E v )a cc (v^E v )dE v (18) 

» Br(LFV)x J &{v e ,E v )<j cc {vn,E v )dE v (19) 

where &lfv is the flux of neutrinos from LFV decays, $(f e ) is the flux of 
electron neutrinos in standard decays and a cc the charged current cross- 
sections. 

For the branching indicated by LSND and 10 19 ^ + decays, we obtain 

N n-,cc ~ 215'000 x Br (LFV) (20) 
= 540 (21) 

for the fiducial detector mass of 10 ton. 
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For muon tracks of a few GeV momentum bent in the magnetic field, we 
expect wrong charge confusion at the level of ~ 10 _3 %. We therefore think that 
the background produced by u„ CC events with mismeasured charge of the muon 
to be less than one event. 

One crucial experimental aspect is the discrimination of muons versus pions. 
In addition, secondary muons from meson decays produce background, espe- 
cially at low energies. In order to assess realistic efficiencies and experimental 
backgrounds, we illustrate results with a detector with the characteristics of the 
ICARUS liquid Argon imaging TPC|37]]. The detector would be on a medium- 
sized LAr vessel surrounded by a dipole magnet. The size of such detector would 
be in the range of 2 x 2 x 6 m 3 for a total mass of about 30 tons. The magnetic 
field would be oriented perpendicular to the drift .E-field and the incoming neu- 
trino direction, in order to bend the charged particles in the direction of the drift 
field, where a resolution in the range of 200 fim is expected from the LAr TPC. 
Given the hadronic interaction length of LAr of A/ = 84 cm, muons which loose 
about 240 MeV/m are distinguished from pions which interact hadronically. 

The effect of simple cuts on background and signal are illustrated in Table ||. 
We assumed in this case a positive muon ring energy and a total of 10 19 standard 
decays. The LFV decay has the branching probability of 2.5 x 10~ 3 , compatible 
with the LSND excess. 

We considered all sources of backgrounds by means of fully generated neutrino 
events. Large event samples were produced for all neutrino species with the 
proper energy distribution with the help of an event generator [38] which includes 
all exclusive final states and a realistic treatment of the low energy region. 

A muon candidate is identified as a track which stops without interacting in 
the Argon or leaves the detector vessel before interacting or stopping. Hence, 
charged and neutral current backgrounds come from hadrons escaping the de- 
tector without interacting or or from actual muons from meson decays. For the 
charged current background, no veto on the positive muon has been included. 

In Table ||, we observe that the presence of LFV decays will produce an 
excess of events with candidates. The visible energy and the candidate muon 
energy spectra are shown in Figure [3|. As expected, the excess is most visible 
at the highest end of the muon energy distribution, since background comes 
from misidentified soft hadrons. A simple cut on the muon candidate energy 
can be used to largely suppress background to negligible level, while keeping a 
large fraction of the LFV signal events, as shown in the Table |2[ The visible 
energy distribution of events for = 2 GeV with a cut on the muon candidate 
momentum of 1.1 GeV (efficiency for LFV of 25%) is shown in Figure ||. The 
statistics is clearly sufficient to constrain any theoretical prediction of the energy 
distribution of the LFV neutrinos. 



7 Results for fi —> e + v e + vt 

The experimental signal consists in the appearance of final state electrons, while 
standard events have positrons in the final state. Therefore, to experimentally 
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Table 3: Effect of cuts on background and signal. We assumed a negative muon ring 
energy E„- of 1, 2 and 5 GeV and a total of 10 19 standard decays. The lepton number 
violating decay has the branching probability of 2.5 x 10~ 3 . E e is the energy of the 
identified electron in the event. 



detect reaction it is mandatory to envisage an experiment endowed with 

charge discrimination capabilities for electrons. This is a true experimental chal- 
lenge given the short radiation lengths in dense targets. 

In order to evaluate realistic efficiencies and backgrounds, we consider an 
experimental setup similar to the one discussed in the previous section, namely, 
an ICARUS liquid Argon imaging TPC with a magnetic field provided by an 
external dipole magnet. 

Table |3| shows the effect of the cuts applied for a normalization of 10 19 muon 
decays. To compute the expected number of signal events we have taken a 
branching probability of 2.5 x 10~ 3 . Given the low muon energies considered, 
most of the events will be quasielastic. We thus require a final state configuration 
containing an electron and a reconstructed proton^] and no additional hadronic 
particles. These criteria reduce the quasielastic background (where we expect 
a neutron rather than a proton in the final state) by almost three orders of 
magnitude, while keeping more than 50% of the signal. The NC background, 
where electron candidates come from 7r° conversion, is in general soft. After 

3 We assume a kinetic threshold of T > 50 MeV in order to detect a proton. 



10 



a cut on the electron candidate momentum, this kind of background becomes 
negligible. 

As mentioned before, charge discrimination between electrons and positrons is 
a must for this kind of search. Preliminary studies show[39] that for a LAr detec- 
tor immerse in a IT magnetic field, a fit to the direction of the electromagnetic 
shower could provide a good determination of the charge of leading electrons. 
When loose criteria are applied, the expected charge contamination from z/ e CC 
amounts to 2% for an electron efficiency of 25%. Tighter requirements reduce 
the charge confusion to the per mil level for an electron identification efficiency 
close to 10%. Table [3| shows that applying loose criteria in the determination of 
the lepton charge is enough to eliminate v e CC background. 



8 Experimental sensitivities 

Table |] shows the expected sensitivities in case a negative result is found as a 
function of the muon energy circulating in the accumulation ring. Three different 
normalizations have been considered: 10 18 ,10 19 and 10 20 muons. For comparison 
we also show the current best limit on the LFV decay /i + — > e + u e u fJ- [[|. 

For a statistics corresponding to 10 19 muon decays we could improve the 
present sensitivity by more than two orders of magnitude. Three orders of mag- 
nitude can be reached for 10 20 muons. Therefore the interpretation of the LSND 
excess in terms of anomalous muon decay that violates lepton flavor and/or total 
lepton number could be thoroughly experimentally tested. 
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Table 4: Achievable limits in case of negative result at the 90%C.L. for /i + — > e^VfV^ 
and pT — > e~v e vn decays with a 10 ton detector for three different number of muon 
decays. 



9 Conclusion 

A negative result from the MiniBOONE experiment would indicate that the 
neutrino flavor oscillation is not the correct hypothesis to explain the excess seen 
in LSND. It would however not contradict other possible non-flavor-oscillation 
interpretations of the effect. 
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In particular, LFV and LVN decays could play a role in the interpretation 
of the LSND excess. A better understanding of these processes would then be 
particularly relevant, if not mandatory. 

A neutrino factory is an ideal machine to probe such anomalous decays of 
the muon. The pure initial state beam allows to look for these decays without 
intrinsic beam contamination. 
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Figure 1: Expected neutrino fluxes for the standard, lepton flavor and lepton number 
violating /i~ decays as a function of the neutrino energy in the center of mass reference 
frame. 
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Figure 2: Neutrino events energy spectrum for muon ring energy of 1, 2 and 5 GeV. 
The Up (line) and v e (dashed) events come from standard n~ decays, while (dotted) 
and v e (dash-dotted) events come from standard /i + decays. 
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Figure 3: Visible energy (upper) and candidate muon (lower) distribution for LFV 
decays (see text) normalized to LSND excess and 10 19 positive muon decays for 10 ton 
detector. The background from neutral current using the characteristics of an ICARUS 
LAr TPC is also shown. 
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Figure 4: Visible energy distribution (see Figure ||) after a cut on the candidate muon 
momentum for LFV signal (see text) and backgrounds expected in a detector with the 
characteristics of an ICARUS LAr TPC. 
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